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1. INTRODUCTION 

In the medical field, the image is nowadays an indispensable support for the diagnosis, the 
establishment of treatment protocols, patient monitoring, and operative preparation. These images must be 
retained for a certain period, the storage required for archiving all data is constantly changing. The images 
are often consulted locally, but they can also be remotely through networks with bandwidth limited. The use 
of compression so quickly revealed to be essential, both to facilitate storage or remote browsing in these 
masses of data [1]. 

The wavelet packet decomposition is an extension of the transforms discrete wavelet to select the 
paving of the time-frequency plane. This choice is made through the selection of a base wavelet packet .In 
general; the base is selected according to the processed signal according to a criterion satisfying the 
constraints application. This database will be called the best basis [2]. Our objective is based on a 
comparative study between the different images types to which it will apply the wavelet packet transform 
with progressive encoder SPHIT. The estimation and judgment of the compressed image quality is given by 
the evaluation parameters (PSNR, MSSIM and FIV). 


2. WAVELET PACKET 

The decomposition of the signal wavelet packet is a generalization of multi-resolution analysis. 
Spaces approximations are always cut in the same way the details of spaces too are direct sum of two sub 
inferior temporal resolutions space (divided by 2). The packets thus wavelet decomposition of a signal is an 
extension of the transforms discrete wavelet [3]. Figure 1 represents the pyramidal algorithm extended for 
obtaining the coefficients. As for the fast calculation of the coefficients of the wavelet transform of a signal, 
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one proceeds by successive filtering and decimation of the signal. Here, the detail coefficients are also 
decomposed. 


Figure 1. Scheme of the decomposition algorithm wavelet packet of a signal 


gı and h, are the quadratic filter mirror QMF filters associated with the scale ø and wavelet y 


functions. C;j are called wavelet packets. Each packet C;j contains n/ j?coefficients in the framework of the 
decomposition of a signal of length n. Coefficients of wavelet packets are rated Cim (K) , where j is the level 
of resolution, m is the spectral band and k is the index of translation. They are obtained by the decomposition 
of the signal on the bases generated by W,, functions [4]: 


Cim(k) = < x(t), 27 2?Wa (27t — k) > 
Wm = 21/2 Uk Zk Win (2t — k) (1) 
Wom+1(t) = 21/2 Yih, Wy (2t — k) 


where, Wo corresponds to the @ function, and W, corresponds to the y function 


2.1. Redundancy of information 
Figure 2 shows the binary tree decomposition of wavelet packets to one of a signal size. 


Figure 2. Binary tree decomposition of a signal into wavelet packets 


At each level of the tree, all the signal information is shown. At a given level of decomposition, 
corresponds to a regular frequency division. At the last level, each frequency is represented by a point, we 
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lost all the temporary information then we have purely frequency decomposition of the signal. Figure 3 
shows the paving of the time-frequency plane for a signal x © R? for each level of the tree [5]. 


LL 


a- Niveau O b- Niveau 1 


c- Niveau 2 d- Niveau 3 


Figure 3. Paving time frequencies corresponding to each level of the decomposition tree 


3. WAVELET PACKET BASE 

Binary tree decomposition wavelet packet thus provides a highly redundant representation of the 
signal. If one wants to work with a non-redundant representation, he has to choose a packet-based that is to 
say a Set of nodes in the tree whose projection in the time-frequency space forms a partition. This is derived 
from the notion of eligible tree. A qualified tree is composed of nodes with 0 or 2 strings; the base will 
consist of all nodes with no string. Considering the frequency division induced decomposition wavelet packet 
this goes back to cover the frequency axis without overlapping. 
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Figure 4. Examples of wavelet packet bases and paving time frequencies corresponding 


Figure 4 shows two examples of wavelet packet bases .The nodes selected in the basic constitution 
are surrounded by squares. We see that if we prune the trees at the selected nodes, eligible trees are obtained. 
The selected packages can form a partition of the frequency axis (horizontal on the tree). The paving of the 
time-frequency induced by the foundation is shown for a signal x © RÈ [6]. 

The packet-based wavelet decomposition from the principle of the AMR, we can reconstruct the 
signal using reconstruction filters associated with decomposition filters. If W,, is the dual basis functions 
associated with Wm, and B is the set of indices {j, m} selected nodes in a packet wavelet basis [4], then 


X= Lgmies Lk Cm K) se Wm (27t— k) (2) 


4. GENERALIZATION TO IMAGES 

The decomposition of an image on a wavelet packet base follows the same principle as the 
decomposition of a one-dimensional signal. We generalize the wavelet decomposition algorithm filtering and 
decimating coefficients details and approximations. In the case of images are act on the lines and columns, it 
will then quadtree decomposition wavelet packets represented by the Figure 5. 

Representation in the Quaternary wavelet packet, each packet C,;j; corresponding to a node of the 
tree contains information regarding the whole image in the frequency band indexed by (i; j) and whose size is 
determined by the resolution level P. The coefficients of the packet are rated Cy;; (k, 1). This decomposition 
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can be interpreted as the decomposition into sub-bands of the image, with increasing frequency 
resolutions [7]. 


Figure 5.Quadtree decomposition wavelet packet 


5. THE ENCODING ALGORITHM SPIHT (SET PARTITIONING IN HIERARCHICAL TREE) 

The SPIHT algorithm [8] involves a list of insignificant sets (LSP), a list of insignificant coefficients 
(LIS) and a list of significant coefficients (LIP). In this work we applied the modified SPIHT algorithm 
which is the most adapted with wavelet packets [9]. The modified SPIHT has the same principle as the 
SPIHT encoder; so for easier representation, the entire sub-band is treated as a node in a wavelet tree. The 
decomposition level of a certain sub-band is the same as that of the dyadic wavelet decomposition subband 
from which it is obtained by further decomposition (i.e. the same level does not imply the same scale - just 
the same position on a dyadic wavelet tree). 

The natural parent (NP) is a subband whose child is on the next finer scale and after lower wavelet 
decomposition tree level. In addition, a NP's child must have the same relative position as NP, if a relative 
subband position is defined for the dyadic subband from which it is obtained by further decomposition. For 
example, each of the children of T3 in Figure 6 has a child who is in the same relative position. And just to 
be mentioned, all parents in dyadic wavelet decompositions are natural. With regular zerotrees, as in SPIHT, 
each parent coefficient has four of its children on the same spatial location and immediately at the next lower 
level of the tree. In wavelet decomposition trees, the parent-child allocation is complicated by the fact that 
there may be different scale differences between the subbands on the adjacent levels [10]. 

Four different cases of inter-level dependencies are organized. All subbands of the next lower level 
are called child subbands here, although when parental conflicts are resolved, the true parent can be found in 
some higher level. 

a. Subband has within the same relative position child subbands which are on a coarser scale, as it occurred 
in the subtree of T2. 

b. Subband has within the same relative position child subband that is on the same scale, as it occurred in 
the subtree of T2. 

c. Natural parent-child relation - subband has within the same relative position child subbandthat is on a 
next finer scale, as it occurred in the subtree of T3. 

d. Subband's relative position overlaps with the relative position of a child subband that is on a scale more 
than one level finer, as it occurred in the subtree of T3. 


Tı 


T2 


Figure 6. SPIHT packet wavelet algorithm 
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6. EVALUATION PARAMETERS OF THE QUALITY 
The commonly parameter used in image compression is the mean square error MSE. This variable 


(error) is defined by the mean square between the pixel (i, j) of the original image 7 (i J ) , and the pixel (i, j) 


of the reconstructed image I(i, j) [11], [12]. 


1 MIN x 2 
MSE =—— >> j)-1. i) 
UAN m=0n=0 (3) 
The peak signal to noise ratio is given by [13], [14]: 
(okay. 
PSNR = 10- logo [aB] 
MSE (4) 


We then evaluate a new paradigm for assessing the quality of medical images, the similarity index compares 
the brightness, contrast, and structure between each pair of vectors, where the index of structural similarity 
(SSIM) between two signals x and y is given by the following expression: 


SSIM (x, y)= I(x, y)-c(x, y): s(x, y) 


Finally, we need a single global measure of overall image quality that is given by the formula [15], [16]: 


mssm (r, î)= = 5 ssm (r, î,) 
j (5) 


M is the total number of local windows in the image. The values MSSIM exhibit greater consistency with the 
visual quality [17]. The Visual Information fidelity parameter (VIF) quantifies Shannon information that is 
shared between the reference and distorted images with respect to the information contained in the reference 
image itself [18]. VIF test is then evaluated as: 


VIF = ———= (6) 


where, I(X;Y / Z) is the conditional mutual information between X and Y, conditioned to Z; s/is a realization 
of S! for a particular image, the index j runs through all sub-bands in the decomposed image [19]. 


7. RESULTS AND ANALYSIS 

In this work, we are interested in the compression of medical image by the transformed wavelet 
packet (PWT) coupled to the coder SPIHT. For this, we have a MRI image of a brain (axial) [20], of size 
512x512 (grayscale) encoded on 8 bpp recorded through an MRI scanner. We vary the bit rate (Rc) of 0.125 
to 2 bpp and calculate the evaluation parameters Simulation results of the applied algorithm is illustrated in 
the Figure 7. From Figure 7 we note that the applied algorithm provides important values of PSNR, MSSIM 
and VIF for a bit rate of 0.5 bpp. 

To better demonstrate the effectiveness and performance of applying PWT in the medical field, we 
will now make a comparison between different types of transformed DWT and PWT coupled with the 
encoder Progressive SPIHT the lifting structure [21], [22]. For this, we chose the same picture MRI of the 
Figure 7 with a level of decomposition N=4 we vary the bit rate from 0.125 to 2 and calculate the evaluation 
parameters. The results are given in Figure 8. 

To illustrate the effectiveness of the applied algorithm we applied different algorithms mentioned 
previously another IRM2 medical image size 512x512 (grayscale) encoded on 8 bpp. We vary the bit rate of 
0.125 to 2 and calculate the assessment parameters. The results obtained are given in Table 1. We distinguish 
that from a bit rate of 0.5 bpp we obtain suitable values of PSNR, MSSIM and FIV. 
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Axial slice of brain Decomposition by the (a) Rc=0.125 bpp; (b) Re=0.25 bpp; 
(Image original) transform wavelet packets MSSIM=0.78; VIF = 0.27; MSSIM=0.88; VIF = 0.42; 
PSNR = 24.81dB PSNR = 30.69dB 


(c) Re=0.5 bpp; (d) Re=0.75 bpp; (e) Rc=1 bpp; MSSIM= 0.99; (£) Rc=2 bpp; MSSIM=1; 
MSSIM = 0.95; VIF = 0.68; MSSIM=0.98; VIF = 0.80; VIF = 0.88; PSNR = 46.63dB VIF = 0.98; PSNR = 55.19dB 
PSNR = 38.24dB PSNR= 43.14dB 


Figure 7. MRI image compression by CDF9/7 (PWT) + SPIHT 
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Figure 8. PSNR, MSSIM, and VIF Variation of MRI image 


Patient - Male 
Age - 52 years 
Exam performed with CTL coil 
Sag FRFSE-XL T2 4mm 
Patient in serious aggravation 
Plaque of demyelization, multiple sclerosis [20] 


Original image Reconstructed image by CDF 9/7 
(PWT) + SPIHT 
PSNR=38.46dB; MSSIM=0.94; 
VIF=0.70 


Figure 9. MRI2 image compression by CDF9/7 (PWT) + SPIHT (Rc = 0.5 bpp) 
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This study was widespread on three medical images (MRI2, TDM and echography) from the 
medical database GE Medical System [23]. The following Figure 9 and Figure 10 show the results obtained 
after the application of PWT + SPIHT algorithm for a rate of 0.5 bpp. As a result, we can state, from 
Figure 9 and Figure 10, that from 0.5 bpp, the applied algorithm gives an acceptable image quality. 


Original image (TDM) Original image (ECHO) Reconstructed TDM image by Reconstructed ECHO image by 
CDF9/7(PWT) +SPIHT CDF9/7(PWT) +SPIHT 
PSNR = 33.20dB, MSSIM = PSNR=34.43dB, 
0.89 ; VIF = 0.55 MSSIM=0.90 ; VIF=0.50 


Figure 10. TDM and ECHO medical images compression by CDF9/7 (PWT) + SPIHT algorithm 
(Re = 0.5 bpp) 


Table 1. PSNR, MSSIM, and VIF Variation of IRM2 medical image 


Bit rate (bpp) 


Filter 0.125 0.25 0.5 0.75 1 1.5 2 
Parameter 
CDF9/7 PSNR (dB) 23.40 30.26 35.60 37.12 40.86 42.34 46.83 
(FB)+ MSSIM 0.50 0.74 0.90 0.92 0.96 0.97 0.99 
SPIHT VIF 0.18 0.39 0.60 0.67 0.78 0.83 0.92 
CDF9/7 PSNR (dB) 24.00 30.66 36.91 40.44 42.81 46.08 48.85 
(lift)+ MSSIM 0.51 0.80 0.92 0.96 0.98 0.99 0.99 
SPIHT VIF 0.18 0.41 0.64 0.76 0.83 0.90 0.94 
GALLS/3 PSNR (dB) 22.84 29.64 35.41 38.62 40.69 44.11 46.93 
(lift)+ MSSIM 0.52 0.76 0.91 0.96 0.97 0.98 0.99 
SPIHT VIF 0.15 0.37 0.60 0.71 0.77 0.86 0.91 
CDF9/7 PSNR (dB) 22.52 27.10 32.63 34.71 36.88 40.79 41.84 
(lift)+ MSSIM 0.45 0.64 0.84 0.90 0.92 0.96 0.97 
EZW VIF 0.14 0.27 0.48 0.57 0.64 0.77 0.80 
GALLS/3 PSNR (dB) 21.58 26.98 30.52 33.58 35.95 37.79 40.35 
(lift)+ MSSIM 0.47 0.66 0.81 0.87 0.92 0.95 0.97 
EZW VIF 0.12 0.27 0.40 0.52 0.62 0.68 0.76 
CDF9/7 PSNR (dB) 29.28 33.44 38.46 41.45 43.45 46.47 49.09 
(PWT)+ MSSIM 0.77 0.87 0.94 0.97 0.98 0.99 0.99 
SPIHT VIF 0.34 0.51 0.70 0.79 0.84 0.90 0.94 
GALLS/3 PSNR (dB) 28.92 32.96 37.02 39.76 41.49 44.52 47.56 
(PWT)+ MSSIM 0.75 0.87 0.94 0.96 0.97 0.99 0.99 
SPIHT VIF 0.34 0.50 0.64 0.74 0.79 0.86 0.92 


8. CONCLUSION 

In this paper, we applied a compression approach on medical images using wavelet packet coupled 
with the progressive coder SPIHT. This approach was tested on medical images (MRI, TDM and ECHO), the 
results are satisfactory. Interesting compression ratio and good quality of the reconstructed images are 
obtained in comparison with other algorithms (DWT, LIFTING). The applied method allows to preserve fine 
structures and produces images with better quality, which is important for medical applications and this is 
proved by the important values of the evaluation parameters. 
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